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ABSTRACT 

We present H- and A'g-band imaging data resolving the gap in the transitional disk around LkCa 15, revealing 
the surrounding nebulosity. We detect sharp elliptical contours delimiting the nebulosity on the inside as well as 
the outside, consistent with the shape, size, ellipticity, and orientation of starlight reflected from the far-side disk 
wall, whereas the near-side wall is shielded from view by the disk's optically thick bulk. We note that forward- 
scattering of starlight on the near-side disk surface could provide an alternate interpretation of the nebulosity. 
In either case, this discovery provides confirmation of the disk geometry that has been proposed to explain the 
spectral energy distributions (SED) of such systems, comprising an optically thick disk with an inner truncation 
radius of ^A-6 AU enclosing a largely evacuated gap. Our data show an offset of the nebulosity contours along 
the major axis, likely corresponding to a physical pericenter offset of the disk gap. This reinforces the leading 
theory that dynamical clearing by at least one orbiting body is the cause of the gap. Based on evolutionary 
models, our high-contrast imagery imposes an upper limit of 21 A/jup on companions at separations outside of 
O'.'l and of 13 A/jup outside of 0'.'2. Thus, we find that a planetary system around LkCa 15 is the most likely 
explanation for the disk architecture. 

Subject headings: stars: pre-main sequence — circumstellar matter — planetary systems — techniques: high 
angular resolution — stars: individual (LkCa 15) 
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1. INTRODUCTION 

The circumstellar disks of gas and dust around newly 
formed stars are believed to be the birthplaces of giant plan- 
ets. In some protoplanetary disks, evidence of gaps or in- 
ner cavities has been revealed through analysis of the infrared 
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spectral energy distribution ( SED, e.g. [C alvet et alj l2002l) or 
interferometry at infr ared (e.g.|R atzka et al.i2007i) or millime- 
ter wavelengths (e.g.l Brown et al. 2009). These objects have 
been termed "transitional disks", since they are thought to 
represent a transitional state of partial disk dissipation be- 
tween the protoplanetary disk stage and the debris disk stage, 
with a variety of mechanisms proposed to account for cen- 
tral cavities: Photoevaporation, disk i nstabilities, and dynam- 
ical clearing (Alexander e t al.l l2006t IChiang & Murrav-Clavl 
]2007; Bryden et al. 1999). 

One such system is LkCa 15 (K5:Ve, V=l 1.91 mag, 
J=9. 42mag, H=8 6mag , K=8.16mag), a 0.97 Mq, 3-5 Myr 
old JSimon et al.l l200(il) , we ak-line T Tauri s t ar with a gas- 
rich m illimeter-bright disk (P ietu et akl l2007t iHenning et alj 
l2010i). Recently, it was identified as a pre-transitional disk 
dEspaillat et alJ2007h , which consists of an inner disk compo- 
nent near the dust sublimation radius (0.12-0. 15 AU), a wide 
gap, and an outer, dis k beyond 46 A U (Espai llat et alj |2008). 
Spitzer IRS spectra bv lSargent et alj (i2009) reveal a relatively 
low level of dust grain growth. With an expected angular scale 
of 0'.'33, the outer wall is accessible to 8-m class telescopes 
equipped with high-contrast imagers. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Subaru/HiCIAO data 

As part of the ongoing high-contrast imaging survey 
SEEDS (Strategic Exploration of Exoplanets and Disks with 
Subaru/HiCIAO; Tamura 2009), we observed L kCa 15 on De- 
cemb er 26, 2009, with the HiCIAO instrument dHodapp et alj 
120081) on the Subaru Telescope with a field of view of 
20"x20" and a plate scale of 9.5 mas/pixel. The images were 
taken in H band ( 1 .6 |am), and the image rotator was operated 
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Fig. 1. — High-contrast imaging of LkCa 15. The scale bars in all panels 
correspond to 1" = 140 AU in projected separation. The greyscale is linear. 
(a) HiCIAO H-hsind image after conservative LOCI ADI. The inner edge of 
the outer disk is clearly visible. The blacked-out area in the center represents 
the range of separations where the field rotation is insufficient for ADI. The 
negative features within (X'S are areas oversubtracted by LOCI as a reaction 
to the nearby positive features. The brightest part of the nebulosity con'e- 
sponds to H = 12.7mag/arcsec^. (b) The signal-to-noise map of the same 
image, calculated in concentric annuli around the star The outer wings of the 
nebulosity are rendered more visible. The stretch is \—2 .5cr, 2.5cr]. (c) In - 
terferometry image of the disk gap at 1.4 mm taken fromPie tu et alj 120061) . 
The cross represents the major and minor axes of the nebulosity, (d) The 
HiCIAO //-band image after reference PSF subtraction. The circular halo 
out to ~0'.'6 is most likely an artifact from imperfect PSF matching, (e) The 
NACO ii's-band image after reference PSF subtraction, for comparison. 



in pupil-tr acking mode to en able angular differential imag- 
ing (ADI. iMarois et al.ll2006l) . The data set comprises 168 
frames, each containing the co-add of 3 unsaturated exposures 
of 1.39 s, for a total integration time of 700.6 s. The image 
resolution provided by the adaptive optics system A0188 was 



close to the diffraction limit, with a full width at half maxi- 
mum (FWHM) of 55 mas at 0'.'6 natural seeing in H. 

The data were corrected for flatfield and field distortion. 
In order to search for point sources such as planets, ADI 
was then performed on the images using the LOCI algorithm 
(Loca lly Optimized Combination of Images, Lafreniere et al] 
l2007h . This form of ADI is the most powerful high-contrast 
imaging method currently available, as evid enced by re- 
cent direct detections of su bstellar companions (iMarois et alj 
lIOOllT^dinann et al.ll2009 ). 

For two frames to be eligible for mutual point-spread func- 
tion (PSF) subtraction at a given working radius, we require 
the differential field rotation arc at that radius to be at least 
3 FWHM = 165 mas, as opposed to the 0.75 FWHM typi- 
cally used for the detection of point-sources. This protects 
extended structures of up to the given azimuthal size scale 
from self-subtraction. We refer to this technique as "conserva- 
tive LOCI". A circularly symmetric circumstellar structure is 
still eliminated under these conditions, but as our data demon- 
strate, a sufficiently elliptical structure can survive by inter- 
secting the concentric annuli around the star at a steep angle. 

An alternate method to reveal circumstellar nebulosity is 
subtraction of the PSF of a reference star, which avoids disk 
self-subtraction at the price of less effective stellar PSF re- 
moval. We observed the star SAO 76662 in pupil-tracking 
mode immediately after LkCa 15 for this purpose. A neutral 
density filter was used with the adaptive optics wave-front 
sensor to achieve the same PSF quality as with the fainter 
LkCa 15. The exposures were centered, flatfielded, distortion- 
corrected, and directly collapsed with a pixel-wise median. 
We subtracted the resulting reference PSF from each of the 
LkCa 15 frames before de-rotation and co-addition to ensure 
optimal matching of the pupil-stabilized stellar PSFs. 

2.2. Archival VLT/NACO data 

For comparison, we retrieved Xg-band coronagraphic 
imaging data of LkCa 15 taken with VLT/NACO in 2007 
from the ESO Science Archive Facility under the program ID 
280.C-5033(A). The observations were made in field-tracking 
mode and thus could not be used for LOCI. However, a refer- 
ence star is included in the data set, allowing PSF subtraction. 
The camera with 13 mas pixel scale and the four-quadrant 
phase mask coronagraph had been used for this run, produc- 
ing high contrast at a small inner working angle. 

3. RESULTS 

3.1. Imaging of the disk gap 

The i^-band HiCIAO images reveal a crescent-shaped neb- 
ulosity around LkCa 15 after both conservative LOCI reduc- 
tion with a frame selection criterion of 3 FWHM and conven- 
tional reference PSF subtraction (Figure [T). With a diameter 
of ~1'.'2 (« 170 AU), it is elongated along the position an- 
gle ~60°, leaving localized traces in the concentric annuli 
around the star in which LOCI operates. The consistency 
of the LOCI and reference subtraction images demonstrates 
that the nebulosity is elliptical enough to survive conservative 
LOCI largely intact. Although some flux is inevitably lost in 
this process, the elimination of the star's PSF in LOCI is supe- 
rior to conventional reference PSF subtraction, revealing the 
sharp inner and outer edges of the structure clearly. The cres- 
cent is also found in the PSF-subtracted A'g-band images from 
VLT/NACO, supporting the interpretation that it represents 
a real astronomical feature. The gap enclosed by the inner 
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Fig. 2. — Contrast (a) and mass (b) of companions detectable at 5 ct around 
LkCa 15, based on the H-band LOCI image after convolution with a circular 
aperture of 4 pixel (38 mas) in diameter. The contrast values are converted 
into co mpanion mass using the COND evolutionary models by Baraffe et^I] 
120031) assuming a distance of 140 pc and an age of 5 Myr (solid curve), 3 Myr 
(dashed curve), and 1.5 Myr (dotted curve), respectively. The vertical dotted 
lines mark the inner working angle of ADI in this data set. The horizontal 
dotted lines guide the eye in (a) and indicate the deuterium burning limit of 
13.7 Af jup in (b). No companion candidates ai'e detected. 



edge matches the gap predicted from the SED dEspaiUat et alj 
120071) as weU as the locaHzed flux deficit seen in 1 .4 mm in- 
terferometry (jKetu et al. 2006, cf. Figure [T]) in terms of size, 
eflipticity, and position angle. 

The negative areas in the inner field represent oversubtrac- 
tion, which is inevitable in LOCI when strong signals are 
present. The algorithm attempts minimize the root of mean 
squares of the residual image, thus the positive mean offset 
of a signal will be lost, given that it cannot be distinguished 
from the positive mean offset of the stellar PSF halo. In con- 
servative LOCI, the frame selection criterion that allows the 
survival of azimuthally extended positive flux concentrations 
also produce extended oversub traction areas. This effect ren- 
ders the flux levels in the resulting image unreliable, but pre- 
serves sharp edges in the original flux distribution. Similarly, 
the location of a point source in regular LOCI is reliable, but 
its photometry is not, requiring the determination of a correc- 
tion factor. 

In its brightest pixel, the crescent in the LOCI image 
reaches a peak intensity of H = 12.7mag/arcsec^, which is 
to be taken as a conservative lower limit. Reference PSF sub- 
traction is expected to retain all of the disk flux, but the cen- 
tral part of the image is dominated by strong residuals from 
the PSF subtraction and therefore does not provide a better 
photometric measurement. 

The small extended signal 0.35" to the southeast of the star 
in the iJ-band images likely corresponds to a real physical 
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Fig . 3 . — Profile of the disk wall edges derived from an 1 1 pixel wide strip 
cut from the LOCI image along the perceived major axis of the ellipses in 
Figure |4] (position angle —32.8°), collapsed along the strip width and then 
median-smoothed on the scale of 5 pixels (~ 1 FWHM), shown as a black 
curve. Median smoothing conserves edges well while removing pixel-to- 
pixel noise. The red vertical curve in the center is the bisector of the two op- 
posed slopes of the profile, demonstrating the asymmetry. The orange curve 
shows the same smoothed profile taken from the LOCI image after reseating 
with a map of the squared distance r^ from the star. The dashed horizontal 
fines guide the eye. 

feature, given that it is visible both after LOCI processing and 
after reference PSF subtraction (Figures[Tt,[T]l). 

3.2. Constraints on point sources 
3.3. Pericenter ojfset 

Our LOCI imaging allows us to set upper limits on the point 
sources present in the vicinity of LkCa 15. Figure |2] shows 
the companion mass we can exclude at the 5 a level as a func- 
tion of separation, assuming the COND evolutionary models 
by Baraffe et al. (2003), a distance of 140 pc, and an age of 3- 
5 Mvr , Simon et aL (.2000,) . Since L kCa 15 is par t of the Tau- 
rus star-forming region, for which IWatson et al.l ( 120091) give 
an age of 1-2 Myr, we also plot a detectable mass curve for 
the age of 1 .5 Myr The image was convolved with a circular 
aperture 4 pixel in diameter, and the noise level is calculated 
as the standard deviation in concentric annuli. We compen- 
sate for the expected flux loss due to partial self-subtraction 
by implanting test point sources in the raw frames and mea- 
suring how they are affected by the LOCI algorithm. We note 
that the bright nebulosity dominates the noise level at all sep- 
arations, thus the calculated upper limit is likely conserva- 
tive, overestimating the true residual speckle noise level. We 
do not detect any statistically significant signals that could be 
c onsidered companion candidates. 

iBonavita et al.l (1201 0) present a more detailed discussion of 
point-source constraints around LkCa 15, including alternate 
evolutionary model assumptions, based on the NACO A's- 
band data. They reach the same contrast levels as we do. 

We find a likely positional offset of the inner and outer 
boundaries of the nebulosity from the star along the system 
semi-major axis (Figures [3] |4^). Furthermore, our inner and 
outer fitted ellipses are rotated by —4° and —3° wit h respect 
to the position angle of 150.7° in lPietuetaP (l2007h . with an 
estimated fitting error of 2°. 

In order to quantify the observed offset, we cut a strip with 
a width of 1 1 pixels = O'.'IO centered on the star and oriented 
along the major axis of our fitted ellipses from the LOCI im- 
age, assuming a —3.5° offset from the 1 .4 mm position angle. 
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Since the disk does not suffer from foreshortening along its 
apparent major axis, the distances along the strip are a direct 
measure of physical distance. We collapse the strip into a one- 
dimensional profile, plotted in Figure |3] 

The intensity drops down to zero over a distance of 54 mas 
(« 8 AU, 1 .0 FWHM) on the West side and 8 1 mas (« 1 1 AU, 
1 .5 FWHM) on the East side of the gap, indicating that the 
edge is sharp at or below the image resolution. The bisecting 
curve between the two opposed gap edges reveals a system- 
atic offset from the star's position by 64 mas (9 AU) with a 
standard deviation of 6 mas (0.9 AU). Since the star is unsat- 
urated, its centering accuracy (^0.2 pixels = 2 mas) does not 
contribute significantly to the uncertainty. The distance of the 
gap edges from the star are 345 mas = 48 AU on the left side 
and 447 mas = 63 AU on the right side at an accuracy of 9 mas 
(1.3 AU), consistent with the 46 AU radius derived from mil- 
limeter interferometry by Pietu et al.l (120061) . 

After scaling each pixel in the strip by r^ before collaps- 
ing the strip, we find that the nebulosity appears at roughly 
the same brightness on both sides of the gap, suggesting a 
r^^ dependency of the unsealed flux levels. This is consistent 
with the assumption that we are looking at reflected light from 
material at varying distances r from the star 

Thus, the inner edge of the outer disk around LkCa 1 5 likely 
features a pericenter offs et comparable to those observed in 
the di sks of HP 14252 7 (iFukagawa et al.ll2006l) and Fomal- 
haut (iKalas et al.ll2005h . 

4. DISCUSSION 

4.1. Forward scattering scenario 

One possible explanation of the observed bright nebulos- 
ity is forward-scattering, which is commonly invoked to ex- 
plain the brightnes s asymmetries in disk s urfaces seen in 
reflec ted light (e.g. [Weinberger et al.l Il999t IFukagawa et al.l 
12006*). Forward-scattering on large dust grains can be sev- 
eral times as efficient as backward-scattering; as a result, the 
near-side surface of such a disk appears brighter than the far 
side. This requires the outer disk surface to be illuminated 
by the star. The firs t of the two SED-compliant models in 
iMulders et al.l (120101) includes a forward-scattering disk sur- 
face; furthermore, the far-side gap wall is shadowed by an 
optically thick inner disk, rendering it hard to detect in di- 
rect imaging (Figure |4};). Note that the size and shape of the 
dust grains i n the wall a r e not well known. This scenario is 
supported bv lPietu et al.l (12007^) . whose orientation and incli- 
nation values for the LkCa 15 disk suggest that the northwest 
side is the near side, as well as the fact that the crescent seen 
after reference PSF subtraction features a bright spot along 
the minor axis. 

4.2. Illuminated wall scenario 

Another explanation for the nebulosity is that it represents 
the inner wall of the outer disk on the far side of the star, illu- 
minated by the star a nd viewed directly through the disk gap 
dEspaillat et al.ll2008i Figure |4]3). The illuminated surface of 
the near-side wall, on the other hand, is blocked from sight by 
the bulk of the optically thick disk. Furthermore, the wall is 
high enough to cast the outer surface of the disk into shadow, 
suppressing a forward-sc attering signa t ure. T his corresponds 
to the second model bv .Mulders et al.l (120101) . The morphol- 
ogy of the nebulosity favors this scenario, given that (1) there 
is a sharp outer edge roughly parallel to the inner edge, (2) 
the nebulosity reaches across the major axis, "embracing" the 




Fig. 4. — (a) Ellipse fits to the inner and outer boundaries of the scattered 
light nebulosity seen in the HiCIAO //-band LOCI image after median fil- 
tering on the spatial scale of 5 pixels fa 1 FWHM and derotation by —29.3° 
(based on the position angle of 150.7° in Pietu et al. 2007). The inner (or- 
ange) and outer (red) ellipses ai'e offset from the stai' along the major axis by 
51 mas and 57 mas and rotated by —4° and —3°, respectively. Their centers 
are marked by orange and red plus signs, respectively, while the star's posi- 
tion is indicated by a white plus sign. The dotted lines delimit the area on 
which the quantitative analysis in Figure[3]is based (offset angle —3.5°). (b) 
Sketch of the illuminated wall scenaiio, taken from Espaillat et al. t2003), 
based solely on the SED and milHmeter interferometry. The inner disk is 
not to scale. The illuminated disk wall on the far side (fight grey) is di- 
rectly visible while the near-side wall blocks its bright side from view, (c) 
_ff-band image of the forward scattering scenario including anisotropic scat- 
tering, from the simulation presented in Mulders et al. (2010). The near-side 
disk surface (top) appears bright due to efficient forward scattering, whereas 
the far-side disk wall (bottom) is mostly shadowed by the inner disk (center, 
not resolved), reducing the wall image to two thin parallel arcs. 

Star, and (3) the nebulosity is wider along the major axis than 
along the minor, all of which are expected for an illuminated 
tapered disk wall, but not for forward scattering on the disk 
surface. Thus, we consider the illuminated wall scenario the 
more likely explanation until further data become available. 

4.3. Gap formation mechanism 

Both scenarios presented above must invoke a disk gap to 
explain the sharp inner edge of the observed nebulosity. Our 
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data therefore prove the vaHdity of the gapped disk model pre- 
dicted from the SED. 

Several mechanisms have been proposed to produce cav- 
ities or gaps in protoplanetar y disks. Magneto-rotational in- 
stability (MRl) as described in lChiang &Murrav-Clavl(l2007h 
can be excluded for LkCa 15, since it acts on all disk compo- 
nents that can be ionized by direct X-ray illumination from 
the star, thus the inner dust component of LkCa 15 known 
from the SED could not have survived. Furthermore, the stel- 
lar mas s of 0.97 Mp^ and stellar accretion rate of 2.4 • 10^^ 
MQ/yr (lEspaillat et a l. 2007) would require a viscosity pa- 
rameter a w 0.0007, one order of magnitude below the range 
considered in the publication. While photoevaporation can in 
principl e evacuate gaps as large as that of LkCa 15 in simu- 
lation s ([Alexander etal.1 120061 iGorti et al.ll2009t lOwen etal] 
1201 Oh . those require that the inner disk drain away entirely 
before the gap can grow beyond 1-10 AU. The presence of an 
inner dust disk around L kCa 15 does not fit this scenario. 

Esp aillat et al.l (l2008h therefore conclude that dynamical 
clearing by one or more orbiting bodies is the most plausi- 
ble cause of the disk gap. Our detection of a likely pericen- 
ter offset in the LkCa 15 disk gap reinforces this argument. 
Dynamical sculpting is commonly invo ked to explain warps 
and eccentricitie s in debris disks (e.g. iRoques et al.l Il994t 
iKalas et al.l2005h . and also produces off-centered gaps in sim- 
ulations of protoplanetary disks (e.g. lMarzari et alj|2010l) . In 
contrast, neither MRI nor photoevaporation have been pre- 
dicted to cause such effects. 

In the case of LkCa 15, close stellar companions are 
excluded by long -baseline interferometry ( Pott et al. 2010) . 
while our data conservatively exclude companions more mas- 
sive than 21Afjup exterior to O'.'l (14AU), and afl but 



planet ary-mass bodies exterior to 0'.'2 (28 AU). iPietu et al] 
( 120061) note that a 5-10 Mjup body orbiting at 30 AU (flie 
equivalent of Neptune's orbit) has a sufficiently large Hill 
sphere to dynamically produce a wall at 50 AU. Less massive 
bodies orbiting closer to the wall can have similar dynamical 
effects. iLubow et al . (1999) note that bodies more massive 
than 6 A/jup will suppress accretion onto a young Solar ana- 
log to levels below that still present in the LkCa 15 system. 
Our data indicate that any bodies in the disk of LkCa 15 at 
r > 0'.'2 (28 AU) from the primary in projection must have 
planetary masses. 

We therefore find dynamical clearing by one or more plan- 
ets to be the most likely cause of LkCa 15's disk gap. Given 
the Sun-like mass of the young star (0.97 M©), we might in 
fact be looking at a Solar System analog in the making. 
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